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摘要 :基于 Gumbel 极 值 1 型 分 布 ,对 在 线 检测 原始 数据 进行 处 理 ,建立 油气 管道 最 大 腐蚀 深度 预测 模型 ,确定 
整 条 油气 腐蚀 管道 可 能 存在 的 最 大 腐蚀 深度 ; 以 此 建立 管道 剩余 壁 厚 模型 ,来 预测 油气 管道 的 剩余 寿命 ;最 后 
以 国内 的 某 一 管道 为 例 ,将 该 模型 应 用 于 整 条 管道 的 剩余 寿命 预测 中 ,预测 结果 为 19.95 a, 表明 该 方法 预测 的 
结果 是 合理 的 。 
关键 词 :油气 管道 腐蚀 最 大 腐蚀 深度 剩余 寿命 预测 Gumbel 极 值 I 型 分 布 Markov 链 Monte Carlo 模拟 
中 图 分 类 号 :TE88 ”文献 标识 码 :A 文章 编号 :1005-4537(2016)04-0370-05 


Residual Life Prediction of Buried Oil and Gas Pipelines 
Based on Gumbel Extreme Value Type I Distribution 


ZHANG Xinsheng, CAO Naining, LI Yayun 
School of Management, Xi'an University of Architecture and Technology, Xi'an 710055, China 


Abstract: Original corrosion data, which were acquired by online detecting a buried oil and gas 
pipeline, were processed in accord with the Gumbel extreme value type I distribution, then a model 
for predicting the maximum corrosion depth is built to determine the possible maximum corrosion 
depth of the whole pipeline with the above method, a model for predicting the residual thickness of 
pipe wall is established to predict the residual life of oil and gas pipelines. In the end, as an exam- 
ple, the model was used to predict the residual life of a domestic pipeline in service, it is ca 19.95 a, 
which indicates that the predict result is correct and reasonable. 
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1 前 言 大 的 经 济 损失 。 因 此 ,研究 腐蚀 对 管道 结构 完整 性 
近年 , 随 着 油气 管道 使 用 年 限 的 增加 , 越 来 越 多 的 危害 程度 , 预测 油气 管道 剩余 寿命 确保 管道 安全 

的 油气 管道 带 缺 陷 工作 ,特别 是 由 于 腐蚀 造成 的 管 。 ”运行 其 有 时 要 的 研究 音义。 
道 缺 陷 越 来 越 严重 , 管道 腐蚀 失效 引起 的 管道 穿孔 、 前 ,国内 外 学 者 对 油气 管道 剩余 寿命 的 预测 
泄漏 事故 频率 也 在 不 断 上 升 ,给 社会 和 企业 造成 巨 。 进行 了 研究 ,如 国外 学 者 Kleiner" 利 用 Markov 模 型 
定稿 日 期 ;2015-12-10 建立 了 管道 寿命 预测 模型 ; 俞 树 荣 等 中 采用 Monte 
基金 项 目 : 国家 自然 科学 基金 项 目 (61271278), 陕西 省 社会 科学 基 Carlo 方 法 计算 以 Shell-92 确定 性 模型 为 基础 所 建立 
金 项 目 (2015R13) 和 陕西 省 自然 科学 基金 项 目 (2016JM6 的 管道 广 蚀 剩余 寿命 预测 概率 模型 ,计算 出 了 管道 
本 。 ee 剩余 寿命 ; 陈 典 斌 等 @ 以 灰色 理论 GM (1,1) 模型 为 
张 新 生 a 基础 ,建立 腐蚀 管道 可 靠 性 模型 并 使 用 VB 进行 腐 
管道 风险 评估 与 管理 及 智能 信息 处 理 等 蚀 管 道 剩余 寿命 分 析 。 采 用 这 些 研究 方法 得 到 的 结 
DOI: 10.11902/1005.4537.2015.152 果 与 实际 值 较 接 近 , 但 是 仍 存在 不 足 。 油 气管 道 的 
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窝 蚀 速率 并 不 是 固定 的 ,而 是 随 着 时 间 变 化 的 ， 
此 , 认为 腐蚀 速率 为 固定 值 预测 出 的 管道 剩余 寿命 
还 是 存在 较 大 误差 的 。 
本 文 将 基于 Gumbel 极 值 I 型 分 布 的 统计 方法 构 
建 油气 管道 剩余 寿命 预测 模型 , 通过 分 析 已 知 管道 
腐蚀 在 线 检测 数据 , 建立 油气 管道 的 最 大 腐蚀 深度 
预测 模型 并 进行 预测 , 最终 预 测 出 油气 管道 剩余 
寿命 。 
2 极 值 I 型 分 布 

对 于 有 腐蚀 缺陷 的 管道 , 缺陷 类 型 通常 为 均 
腐蚀 、 局 部 腐蚀 和 点 刨 甸 。 而 影响 管道 使 用 寿命 
主要 因素 是 点 腐蚀 , 因为 引起 管道 穿孔 破坏 的 是 最 
深 的 蚀 孔 ,一 旦 管道 穿孔 , 管道 就 发 生 泄漏 失效 事 
故 ,因此 ,最 大 的 腐蚀 坑 深 度 是 评定 穿孔 破坏 的 重要 
指标 外 。 由 于 油气 管道 工作 环境 复杂 , 其 腐蚀 过 程 
的 规律 难以 掌握 ,通过 建立 统一 的 数学 公式 来 表示 
各 腐蚀 因素 的 影响 规律 是 比较 困难 的 路。 统计 方法 
是 通过 统计 分 析 管 道 最 大 腐蚀 深度 , 通过 计算 腐蚀 
失效 概率 来 预测 管道 剩余 寿命 后。 极 值 I 型 分 布 是 
观测 值 中 最 大 值 或 最 小 值 的 概率 分 布 ,是 概率 统计 
方法 中 的 一 个 重要 分 支 ,广泛 应 用 于 许多 工程 实际 
问题 的 分 析 、 评 佑 中 。 

对 于 观测 值 的 最 大 值 或 最 小 值 的 集合 进行 统计 
的 极 值 1 型 分 布 形式 如 下 : 


Fi(y) = exp -emp 人 - 2 2)) XE (-%,%); 
Ae(-%,%);a>0 (1) 


此 即 Gumbel 极 值 I 型 分 布 中 的 最 大 值 分 布 , 又 称 为 
Gumbel 极 值 I 型 分 布 ,其 中 , Fi(W 为 最 大 值 不 超过 x 
的 概率 ,x 为 某 一 随机 变量 的 最 大 值 ; a 为 尺度 参数 ， 
4 为 位 置 参数 。 


Fi(y) = 工 一 ep| -emp(- 此 2 )} XE (-%,%); 
Ae(-%,%);a>0 (2) 


此 即 Gumbel 极 值 I 型 分 布 中 的 最 小 值 分 布 , 其 中 
五 (7 为 最 小 值 不 超过 x 的 概率 ,x 为 某 一 随机 变量 的 
最 小 值 。 

管道 腐蚀 缺陷 的 概率 分 布 有 一 定 的 规律 ,例如 ， 
应 力 腐蚀 破裂 寿命 符合 Weibull 分 布 , 管道 腐蚀 增长 
量 符合 逆 Gaussian 分 布 ,管道 腐蚀 点 数 符合 Possion 
分 布 等 中 。 学 者 经 过 大 量 的 实验 数据 分 析 中 认为 ， 
油气 管道 的 局 部 腐蚀 的 最 大 腐蚀 点 的 概率 分 布 服 从 
Gumbel 极 值 I 型 分 布 函数 。 
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3 基于 极 值 I 型 分 布 的 管道 最 大 腐蚀 深度 预测 
3.1 极 值 I 型 分 布 应 用 于 埋 地 油气 管道 腐蚀 剩余 寿命 
的 预测 

极 值 I 型 分 布 用 于 油气 管道 腐蚀 剩余 寿命 预测 
的 基本 步骤 如 下 :() 获取 腐蚀 严重 区 段 的 管道 腐蚀 
数据 ;(2) 将 获得 的 数据 通过 数理 统计 方法 计算 最 大 
腐蚀 深度 的 分 布 概率 利用 图 形 检验 的 方法 来 检验 最 
大 腐蚀 深度 符合 极 值 分 布 是 否 合理 ;(3) 若 符合 极 值 
分 布 , 则 构造 统计 分 析 模 型 ; (4) 利用 Markov 链 
Monte Carlo (MCMO) 仿真 估计 模型 参数 ; (5) 基于 
极 值 I 型 分 布 获得 整 条 管道 可 能 的 最 大 腐蚀 深度 ; 
(6) 由 APIRP5795 推 荐 的 公式 求 出 腐蚀 管道 的 临界 
腐蚀 深度 ; (7) 最 终 基 于 所 得 最 大 腐蚀 深度 和 临界 腐 
刨 深度 建立 剩余 寿命 预测 模型 来 预测 油气 管道 的 剩 
余 寿 命 。 油 气管 道 腐蚀 剩余 寿命 预测 的 基本 步骤 如 
图 1 所 示 。 
3.2 油气 管道 最 大 腐蚀 深度 预测 模型 

管道 的 局 部 腐蚀 最 大 值 的 概率 分 布 符合 Gum- 
bel 极 值 1 型 分 布 , 其 分 布 函数 依据 式 (1) 为 : 


Fi(7)= -ee = i 


(3) 
式 中 ,A() 为 最 大 腐蚀 深度 不 超过 x 的 概率 ;x 为 最 
大 腐蚀 深度 的 随机 变量 ; 4 为 统计 参量 , 表示 概率 密 
度 最 大 的 腐蚀 点 深度 ; a 为 统计 变量 ,表示 腐蚀 点 深 
度 的 平均 值 ; 变 量 代 换 为 产 CCA)/a。 

由 于 油气 管道 埋 在 地 下 ,在 收集 管道 腐蚀 坑 深 
度数 据 时 ,不 可 能 将 整 条 管道 开 挖 ,采集 所 有 腐蚀 信 
息 , 这 样 工程 量 巨大 难以 实现 且 成 本 高 。 通 常情 况 ， 
工程 上 将 管道 外 防腐 层 检测 严重 的 区 域 进行 开 挖 ， 
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图 1 剩余 寿命 预测 步骤 流程 图 


Fig.1 Flow chart of residual life prediction 
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进行 取样 收集 数据 。 这 样 获得 的 最 深 的 腐蚀 深度 存 
在 很 大 的 随机 性 , 能够 反映 整个 管道 的 腐蚀 情况 。 
将 检测 的 数据 进行 统计 处 理 , 得 到 统计 参量 wc 和 1, 然 
后 根据 式 (3) 计算 最 深 的 腐蚀 深度 不 超过 某 一 数值 
的 概率 。 

假设 有 N 个 腐蚀 深度 x 这 些 腐蚀 深度 均 是 来 
自 于 处 于 相同 环境 的 同一 条 管道 的 管 段 上 的 最 大 腐 
刨 深度 点 ,将 和 N 个 在 同样 条 件 下 测 得 的 最 大 腐蚀 坑 
深度 数据 由 小 到 大 排序 , 标 出 序号 7 和 对 应 的 腐蚀 坑 
深度 值 xs, 计算 累积 分 布 函数 : 
F(y,)=i/N+1,(=1,2,3,.…,N) (4) 


将 式 (3) 两 端 同时 取 2 次 对 数 , 得 : 
y= -mm 站 =(%— N/a (5) 
获得 统计 参量 ga 和 4, 代入 到 式 (3) 就 可 以 求 出 


最 大 腐蚀 深度 不 超过 茶 一 数值 的 概率 FX<x)。 以 % 


为 模 些 标 ，-In|-In 志 二 | 为 纵 坐标 ,将 各 点 (x ， 


-| 二 中 画 在 直角 坐标 系 中 ,如 果 这 些 点 呈现 
I i 


线性 关系 , 说 明 这 些 腐蚀 数据 服从 Gumbel 极 值 I 型 
分 布 , 然 后 进行 统计 参量 估计 。 
回归 周期 7(W) 与 累计 概率 (7) 的 关系 为 : 


TW=1/(1-F(»)) (6) 


回归 周期 是 某 测 定 值 x 为 观测 顺序 号 的 期 待 
值 。 即 抽取 入 个 样本 ,其 最 大 值 意 味 着 从 少量 小 段 
的 试 样 最 大 坑 深 推算 出 的 整个 局 部 段 的 坑 深 。 
TW)=L/N (7) 

根据 式 (7) 就 可 以 得 到 整 条 工 长 管道 大 致 的 最 
大 腐蚀 坑 深 度 。 将 式 (6) 求 得 的 RD) 值 代入 式 (5) 
可 求 得 y 值 ;最 后 将 y 值 代入 拟 合 直线 方程 可 以 求 出 
最 大 局 部 腐蚀 坑 深 度 。 管 道 腐蚀 段 工 中 的 最 大 腐蚀 
点 深度 的 预测 值 x, 为 : 

X =X%=Atay (8) 
3.3 预测 模型 参数 的 估计 

利用 MCMC 模拟 对 油气 管道 最 大 腐蚀 深度 预 
测 模 型 的 参数 进行 估计 ,在 无 信息 先 验 分 布下 ,用 
MCMC 专用 软件 包 OpenBUGS 对 预测 模型 参数 进 
行 估计 。 

经 过 对 大 量 的 实验 数据 分 析 研 究 发 现 * 蚊 , 油 
气管 道 的 局 部 腐蚀 的 最 大 腐蚀 点 服从 Gumbel 极 值 I 
型 分 布 函数 , 即 六 Gumbel(w, 办 。 考 虑 到 Gamma 分 
布 能 够 保证 a 和 4 为 正 值 且 可 以 方便 的 作为 一 个 无 信 
息 分 布 , 本 文选 取 Gamma 分 布 作为 a 和 4 的 先 验 分 


布 , 妈 a~Gamma(pi; 91); 4~Gamma(py, 0), 其 中 ,pi(q') 
和 py(q;) 分 别 表 示 a 和 4 的 Gamma 分 布 形状 参数 。 

MCMC 方 法 模拟 的 样本 为 10000 个 ,对 Gumbel 
极 值 1 型 分 布 进行 Gibbs 抽样 ,获得 参数 的 Gibbs 迭 
代 轨 和 迹 (trace) 和 迭代 历史 (history), 根据 大 量 的 参 
数 估计 值 的 趋势 收敛 情况 ,确定 估计 参数 的 准确 性 。 
3.4 腐蚀 管道 的 临界 腐蚀 深度 的 计算 

腐蚀 管道 的 临界 腐蚀 深度 的 计算 可 由 
APIRP579094 推 荐 的 公式 求 出 : 


到 


_ PD 
hs 本 2F0o,E (9) 
Ah=h—h,,. (10) 


式 中 , jw 是 管道 临界 腐蚀 深度 , mm; P 是 管道 运行 
内 压 , MPa; D 是 管道 公称 直径 , mm; 为 设计 系数 ， 
与 管线 所 处 地 区 类 别 相 关 ; os 是 管材 最 小 届 服 强 
度 , MPa; E 是 管材 弹性 模 量 , MPa; 是 管道 公称 壁 
厚 ,mm; Ah 肛 管道 人 允许 腐蚀 深度 , mm 。 
4 油气 管道 腐蚀 剩余 寿命 预测 

实际 油气 管道 的 腐蚀 并 不 是 按照 固定 的 腐蚀 
速度 进行 的 ,一 般 腐 蚀 程 度 会 随 着 时 间 的 增加 而 递 
增 , 但 是 腐蚀 速率 会 随 着 腐蚀 深度 的 增长 而 减 小 。 
根据 腐蚀 深度 与 使 用 年 限 的 关系 ,利用 工程 上 的 经 
验 公式 来 近似 模拟 管道 最 大 腐蚀 深度 与 使 用 年 限 
的 关系 忠 , 预测 管道 的 使 用 寿命 ,从 而 确定 管道 的 剩 
余 寿 命 ,如 下 式 所 示 : 


‘|[ 妾 ] GD) 


由 

式 中 ,# 为 预测 的 管道 设计 使 用 寿命 ,a; 1 为 管道 服役 
时 间 ,a;p 为 t 时 刻 测 到 的 最 大 腐蚀 深度 , mm; c 为 时 
间 常 数 ,由 腐蚀 体系 等 因素 决定 , 只 有 一 个 时 间 测 量 
值 时 ,通常 采用 比较 保守 的 值 为 0.505。 
5 应 用 实例 

辽宁 某 陆 上 油田 集 油 管道 已 经 运行 了 约 10 a， 
管道 规格 为 D273 mmx8(9) mm 油气 管道 材质 为 
X52 钢 , 管道 外 径 D=273 mm, 管道 壁 厚 取 较 小 值 
8 mm, 屈服 强度 wm =290 MPa, 该 管道 工作 压力 产 
0.85 MPa, 不 考虑 地 区 类 别 , 设计 系数 取 1, 即 按照 
100% 届 服 强度 计算 ,管材 安全 系数 取 1.1。 该 管道 
总 长 为 3.6 km。 工 作者 针对 外 防腐 层 检测 腐蚀 严重 
区 域 开 挖 ,分 段 获取 数据 ,检测 时 每 km 作为 一 个 检 
测 段 ,划分 为 21 个 段 ,每 段 随机 检测 60 个 点 ,在 每 个 
检测 段 上 找 出 最 深 的 腐蚀 深度 见 表 1。 
5.1 油气 管道 最 大 腐蚀 深度 模型 的 验证 

取 到 各 段 检 测 的 最 大 腐蚀 深度 代入 式 (5)~(7) 
计算 出 各 点 的 yy 以 x 为 横 坐 标 ,yj 为 纵 坐 标 ,绘制 直 
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角 坐 标 图 并 进行 线性 拟 合 ,如 图 2 所 示 。 线 性 拟 合 
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度 R=0.982 (R 越 接近 1, 说 明 越 接近 线性 ), 说 明 局 
部 腐蚀 管道 中 所 获取 的 最 大 腐蚀 深度 值 与 y 有 线性 
关系 , 这 也 正 说 明了 所 eo 
值 1 型 分 布 。 

5.2 油气 管道 剩余 寿命 预测 

根据 实例 中 提供 的 数据 ,采用 MCMC 方 法 来 估 
计 油 气管 道 最 大 腐蚀 深度 预测 模型 的 参数 , MCMC 
估计 参数 结果 oa=1.323, 三 0.41。 

MCMC 方 法 对 预测 模型 进行 Gibbs 抽样 , 得 到 
预测 模型 参数 的 Gibbs 迭代 轨迹 和 和 办 代 历史 , 其中， 
迭代 轨迹 随机 选择 第 9800~10000 个 样本 进行 迭代 ， 
如 图 3~5 所 示 。 
昌 图 3 和 4 所 示 获 得 的 参数 w 和 41 的 迭代 轨迹 和 
友 代 历史 基本 趋 于 稳定 , 认为 迭代 过 程 已 经 收敛 。 
迭代 收敛 后 ,通过 MCMC 方法 得 到 参数 wc 和 4 的 置 
度 图 形 见 图 5。 两 个 参数 的 分 布 区 间 与 从 代 历史 区 
间 基 本 符合 , 由 此 对 两 参数 做 出 综合 推断 , 从 而 估计 
出 模型 参数 。 

MCMC 方 法 得 到 的 模型 参数 ga 和 4 的 标准 差分 
别 为 0.0977 和 0.0781, 且 和 迭代 轨迹 和 和 迭代 历史 基本 
趋 于 稳定 ,可 知 MCMSC 估计 参数 是 比较 精确 的 。 


表 1 局 部 管道 最 大 腐蚀 深度 


Table 1 Maximum partial corrosion depth 


EH 


QI 


Intervals of maximum Number of maximum 


partial corrosion depth partial corrosion depth /mm 


(0.5, 1.0] 2 
(1.0, 1.5] 10 
(1.5, 2.0] 5 
(2.0, 2.5] 3 
(2.5, 3.0] 1 


yimm 


0.5 1.0 1.5 2.0 2.5 3.0 
Maximum corrosion depth / mm 
图 2 最 大 腐蚀 深度 概率 分 布 
Fig.2 Probability distribution of the maximum corrosion 
depth 


此 ,通过 软件 模拟 ,估计 参数 a=1.323,4=0.41, 则 由 
式 (8) 得 管道 最 大 腐蚀 深度 =1.323+0.41y。 

由 式 (8)~(10) 得 =4.32 mm; 由 式 (9) 得 管道 
临界 腐蚀 深度 和 as=0.44 mm, 由 式 (10) 得 最 大 允许 
腐蚀 深度 Ah =7.56 mm 将 预测 的 最 大 腐蚀 深度 x 
和 Ah 带 入 式 (11) 得 管道 的 可 使 用 寿命 f=30.63 ai; 
整 条 管道 的 腐蚀 剩余 寿命 六 30.63-10.68=19.95 a。 


[ Ge) 


9800 9850 9900 9950 10000 
Sample size 


9800 9850 9900 9950 10000 
Sample size 


图 3 参数 ga 和 4 的 迭代 轨迹 


Fig.3 Iterative trace of parameters a (a) and 4 (b) 


“0 2000 4000 6000 8000 10000 
Sample size 


0 2000 4000 6000 8000 10000 
Sample size 


图 4 参数 a 和 4 的 迭代 历史 
Fig.4 Iterative history of parameters & (a) and4(b) 
0.8 


0.6 


0.4 


Corrosion level 


0.2 


0.75 1.00 1.25 1.50 
Estimation parameter 


(b) 


Corrosion level 


0.25 0.50 0.75 1.00 
Estimation parameter 


图 $ 参数 a 和 4 的 置信 和 度 
Fig.S Confidencial level of parameters w (a) and 4(b) 
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8.0 


t=(7.56/0°t 


6.0 (30.63, 7.56) 


A (10.68, 4.32) 


2.0 


0.0 n 1 1 4 4 1 上 
-1.0 30 70 11.0 15.0 19.0 23.0 27.0 31.0 35.0 


Corrosion time /a 
图 6 管道 剩余 寿命 预测 
Fig.6 Prediction of pipeline residual life 


国外 油气 管道 的 设计 寿命 一 般 为 33 a", 显然 
以 Gumbal 极 信 I 型 分 布 预测 出 最 大 席 蚀 深 度 为 基础 
所 计算 出 的 管道 使 用 时 间 30.63 a 是 合理 的 。 说 明 
通过 利用 概率 统计 的 方法 预测 出 油气 管道 的 剩余 寿 
命 为 19.95 a, 与 该 油气 管道 的 整个 寿命 期 限 比较 , 是 
合理 的 。 


Maximum corrosion depth / mm 


6 结论 
(1) 依据 管道 局 部 挖 坑 在 线 检测 数据 所 得 的 最 


大 腐蚀 深度 概率 分 布 ,与 Gumbel 分 布 能 够 较 好 的 拟 
合 , 油气 管道 的 局 部 腐蚀 的 最 大 腐蚀 深度 的 概率 分 
布 服 从 Gumbel 极 值 I 型 分 布 函数 。 
(2) 用 Gumbel 极 值 I 型 分 布 预测 出 埋 地 油气 
道 的 最 大 腐蚀 深度 , 并 估算 出 管道 的 剩余 寿命 。 
例 表 明 获 取 到 的 预测 值 确定 的 管道 剩余 寿命 是 
理 的 。 
(3) Gumbel 极 值 I 型 分 布 管道 剩余 寿命 预测 模 
型 能 够 较 精 确 的 预测 油气 管道 剩余 寿命 。 
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